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ABSTRACT 

We present observations at 1.2 mm with MAMBO-II of a sample of z <; 2 radio-intermediate obscured 
quasars, as well as CO observations of two sources with the Plateau de Bure Interferometer. The 
typical rms noise achieved by the MAMBO observations is 0.55 mJy beam~^ and 5 out of 21 sources 
(24%) are detected at a significance of > 3a. Stacking all sources leads to a statistical detection of 
(<S'i.2 mm) = 0.96 ± 0.11 mJy and stacking only the non-detections also yields a statistical detection, 
with (51.2 mm) = 0.51 ± 0.13 mJy. At the typical redshift of the sample, z = 2, 1 mJy corresponds 
to a far-infrared luminosity Lfir'~-" 4 x 10^^ Lq. If the far-infrared luminosity is powered entirely 
by star-formation, and not by AGN-heated dust, then the characteristic inferred star-formation rate 
is '^700 M0 yr^^. This far- infrared luminosity implies a dust mass of 3 x 10^ M©, which is 

expected to be distributed on ^kpc scales. We estimate that such large dust masses on kpc scales 
can plausibly cause the obscuration of the quasars. Combining our observations at 1.2 mm with mid- 
and far-infrared data, and additional observations for two objects at 350 fim using SHARC-II, we 
present dust SEDs for our sample and derive a mean SED for our sample. This mean SED is not well 
fitted by clumpy torus models, unless additional extinction and far- infrared re-emission due to cool 
dust are included. This additional extinction can be consistently achieved by the mass of cool dust 
responsible for the far-infrared emission, provided the bulk of the dust is within a radius ~2-3 kpc. 
Comparison of our sample to other samples of z ^ 2 quasars suggests that obscured quasars have, 
on average, higher far-infrared luminosities than unobscured quasars. There is a hint that the host 
galaxies of obscured quasars must have higher cool-dust masses and are therefore often found at an 
earlier evolutionary phase than those of unobscured quasars. For one source at z — 2.767, we detect 
the CO(3-2) transition, with 500^*^ =630±50 mJy km s~^, corresponding to Lcnr.s-a^ =3.2x10^ Lr. 



CO(3-2) =O.ZAlU IJQ, 

or a brightness-temperature luminosity of ico(3-2) ^ ^-^ ^ ^^^'^ ^ P*"^- another source at 

z — 4.17, the lack of detection of the CO(4-3) line suggests the line to have a brightness-temperature 
luminosity ^^0(4-3) ^ 1 ^ ^0^" ^ • Under the assumption that in these objects the high-J 

transitions are thermalised, we can estimate the molecular gas contents to be Mb_2 — 1-9 x 10^" Moand 
< 8 X 10^ Mq, respectively. The estimated gas depletion timescales are Tg — 4 Myr and <16 Myr, 
and low gas-to-dust mass ratios of Mg/M^— 19 and < 8 are inferred. These values are at the low end 
but consistent with those of other high-redshift galaxies. 

Subject headings: galaxiesmuclei - galaxies: active - quasars:general - infrared: galaxies - galax- 
ies:starbursts - ISM:molecules 
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1. INTRODUCTION 

Quasars are believed to be powered by supermassive 
black holes (SMBHs) rapidly growing by accreting mat- 
ter at a large fraction of their Eddington rate^*'. The 
accreted matter forms a disk, is heated to temperatures 
of ^ 10^ — 10^ K, and emits thermal emission at optical 
and ultraviolet wavelengths. Clouds of gas ionised by this 
radiation will emit line radiation, and lines which are rel- 
atively close to the black hole have rapid orbits which will 
show broad velocity dispersions (>2000 km s^^), while 
clouds further away will show narrower lines. 

Dust can survive when the equilibrium temperature 
with the ultraviolet photons is ^2000 K. This dust ab- 
sorbs optical and ultraviolet radiation and re-emits it at 
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In this article we use the term quasar for the most powerful 
active galactic nuclei (AGNs). In the case of unobscured AGNs, the 
dividing line is taken at a B-band magnitude of Mb = —23.5 which, 
assuming the typical quasar SED of Elvis et al. (1994), corresponds 
to Lbol 2 X 10^^ Lq. We will therefore also refer to obscured 
AGNs as obscured quasars if they have Lbol ~ 2 X 10^^ Lq. 
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infrared wavelengths characteristic of the dust temper- 
ature. If the line of sight of an observer to the central 
region is blocked by dust, the optical and ultraviolet con- 
tinuum, as well as the broad lines, will not be observable. 
Such a source is known as an obscured quasar. 

While radio-loud obscured quasars have been known 
for a long time, in the form of high-excitation narrow- 
line radio galaxies (see Urry & Padovani 1995, for a 
review), the population of radio-quiet obscured quasars 
had remained elusive, with only a few individual objects 
known. This has changed recently, when large numbers 
of these sources have been found in the spectroscopic 
database of the Sloan Digital Sky Survey (Zakamska 
et al. 2003). Samples of spectroscopically-confirmed ob- 
scured quasars have now also been identified in surveys 
using X-ray selection (e.g. Szokoly et al. 2004; Barger 
et al. 2005) as well as mid-infrared, sometimes combined 
with radio, selection (e.g. Houck et al. 2005; Martmez- 
Sansigre et al. 2005; Lacy et al. 2007). It seems that the 
obscured quasars are at least as common as the the un- 
obscured quasars (e.g. Reyes et al. 2008), and probably 
outnumber these bv a ~2:1 ratio (e.g. Martmez-Sansigre 
et al. 2005, 2008; Lacy et al. 2007; PoUetta et al. 2008). 

There is a large body of evidence supporting obscu- 
ration of quasars as an orientation effect, where ob- 
scured and unobscured quasars are intrinsically identical 
sources, only dust around the accretion disk covers cer- 
tain lines of sight so that the broad emission lines and 
the optical and idtraviolet continuum are not visible (the 
dusty torus of the unified schemes, e.g. Urry & Padovani 
1995). However, dust distributed within the host galaxy 
can also cause obscuration, particularly if the galaxy is at 
an early evolutionary phase and is rich in gas and dust 
(see e.g. Lawrence & Elvis 1982; Sanders et al. 1988; 
Fabian 1999). 

The importance of galaxy-scale dust is evident in re- 
cent studies of quasars where obscuration cannot be 
solely assigned to dust from a torus around the accretion 
disk: in some objects even the narrow-lines are obscured 
(Martmez-Sansigre et al. 2005, 2006b; Rigby et al. 2006). 
The jets emanating from some obscured quasars are face- 
on (Martmez-Sansigre et al. 2006a; Sajina et al. 2007; 
Klockner et al. 2009), something that is not expected for 
quasars obscured by the torus, but can be explained if 
the obscuring dust is on a larger (~kpc) scale. The deep 
silicate absorption features and spectral energy distribu- 
tions (SEDs) suggest foreground extinction is sometimes 
present (Levenson et al. 2007; PoUetta et al. 2008). This 
all provides indirect evidence suggesting obscuring dust 
on ~ kpc scales, something which has been spectroscop- 
ically confirmed through the measurement of the Ha to 
H/3 Balmer decrement (Brand et al. 2007). 

Dust present on kpc scales is expected to be rela- 
tively cool [T ^50 K) and will emit thermally in the 
far- infrared (A ^40 fim), regardless of whether it is 
heated by young stars or by the central active galactic nu- 
cleus (AGN, see e.g. Barvainis 1987; Sanders et al. 1989; 
Siebenmorgen et al. 2004, and Section 3.5). Observing 
this dust around the peak of its emission (aroimd rest- 
frame ^65 /xm) is not optimal, since warmer dust can 
still contribute significantly there. The Rayleigh- Jeans 
tail of the emission, however, will have a much smaller 
contribution from warm dust, and has the additional ad- 
vantage of being observable from the ground, for example 



in the atmospheric windows at 850 fiia or 1.2 mm. 

In this paper we present continuum observations of a 
sample of 21 high-redshift {z ^2) radio-intermediate ob- 
scured quasars at 1.2 mm, combined with other infrared 
and submillimetre data, as well as a search for CO in two 
of the sources. The sample was selected to approximately 
match the break in the unobscured luminosity function 
at z ~ 2 (Mb ~-25.7, Groom et al. 2004), so in terms 
of radiation these quasars represent the energetically- 
dominant population around the peak of the quasar ac- 
tivity. The sources are, however, radio-intermediate in 
that their radio luminosities arc slightly higher than 
those expected from radio-quiet quasars (with Li.4ghz 
~ 10^'' W Hz-i sr-i Martmez-Sansigre et al. 2006a). 
Nevertheless, observations with very large baseline inter- 
ferometry suggest the relative importance of the compact 
cores makes them more similar to the genuinely radio- 
quiet population than to the radio galaxies and radio- 
loud quasars (see Klockner et al. 2009). Section 2 sum- 
marises the observations and data reduction. In Section 3 
we show the resulting fluxes, inferred luminosities, star- 
formation rates, dust masses, characteristic scale of the 
cool dust and discuss the possible obscuration by kpc- 
scale dust. Section 4 shows the results of broad-band 
SED fitting and Section 5 presents the mean SED for our 
sample. In Section 6 we discuss the molecular gas obser- 
vations. We compare our sample to other millimetre or 
submillimetre observations of z ~ 2 quasars in Section 7. 
The discussion and summary are found in Section 8. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. MAMBO continuum observations 

The sample of obscured quasars was observed at a 
wavelength of 1.2 mm (250 GHz) using the 117-element 
Max-Planck Millimetre Bolometer Array (MAMBO-II 
Kreysa et al. 1998, hereafter simply MAMBO), at the 
IRAM Pico dc Velcta 30 m telescope. The observa- 
tions were carried out in "ON-OFF" mode during the 
MAMBO pool at the IRAM 30 m telescope throughout 
various semesters. AMS16 was observed in winter 2005- 
2006, while the remaining objects were observed between 
summer 2007 (Sum07) and winter 2007-2008 (Win07). 
Due to a warning issued by IRAM (see below), the ob- 
jects detected during the Win07 pool were re-observed in 
summer 2008 (Sum08) and winter 2008-2009 (Win08). 

To minimise the effect of atmospheric absorption, the 
science targets were observed whenever possible at ele- 
vations >45 deg, and often at around 60 deg. Observing 
at elevations higher than 75 deg is, however, not possible 
with the 30 m telescope. 

For the first-order removal of the sky, the observations 
make use of a secondary mirror to chop (or "wobble"), 
with a typical azimuthal throw of 32 to 45 arcsec and 
a chopping frequency of 2 Hz. The ON-OFF observa- 
tions are done by nodding, so that the sky-only position 
sampled by the chopping is alternatively on the right or 
on the left of the sky-plus-source position. Each nod 
position (ON or OFF) is observed for one minute, and 
to minimise overheads an ON-OFF-OFF-ON sequence is 
used. Hence integration blocks ("scans") are multiples 
of 4 minutes. 

The ON-OFF observations were carried out centred 
on the most sensitive pixel, number 20, with the atmo- 
spheric opacity T230 GHz ^0.3 and medium-to-low sky 
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noise (^200 mJy bcam^^). Each source was typically 
observed in several blocks of 20 minute integrations. 

The "sky-only" position can sometimes fall on a real as- 
tronomical source. Given that high-redshift AGNs are of- 
ten found to have (sub)millimetre-bright galaxies nearby 
(Stevens et al. 2003), there is a real risk of ruining the 
sky subtraction by chopping onto a bright source. To 
minimise the risk of this happening (and the impact on 
the final data), the chop angle ("wobbler throw") was 
alternatively varied between 32, 35 and 42 arcsec, and 
observations were carried out on different dates and at 
different times of the day. This means the chop in az- 
imuth corresponds to different physical offsets in RA and 
Dec, so the chance of repeatedly falling on the same as- 
tronomical source is minimised. 

Throughout the observing runs, gain calibration was 
performed by observing Neptune, Uranus or Mars, and 
the flux calibration was always found to be accurate to 
<20%. The gain calibration was also monitored regularly 
by observing other sources with known bright millime- 
tre fluxes (>5 Jy). Total power measurements at differ- 
ent elevations were made to infer the atmospheric opac- 
ity ( "skydip" ) . These were performed typically every 2 
hours under stable weather conditions, and at shorter in- 
tervals for less stable conditions. The skydips were gen- 
erally done with the telescope azimuthal angle matched 
to that of the science targets. The focus of the telescope 
was also regularly monitored on bright sources. Before 
the first science target and every 40 minutes the accuracy 
of the telescope pointing was checked using the nearby 
source J1638-I-573 (0.6 Jy). Observations of J1638-F573 
and other bright sources were also used to check for 
anomalous refraction, in which case no science observa- 
tions were made. 

During the observations, the sensitivity of MAMBO in 
a 1 second integration was typically between 35 and 45 
mJy beam""'^. Including typically 60% overheads, the ob- 
servations of this sample of 21 obscured quasars totalled 
approximately 65 hours. 

The data were reduced using the MOPSIC pipeline 
(developed by R. Zylka). MOPSIC removes the atmo- 
spheric emission using the two chopper positions'^ The 
sky noise is measured by taking the weighted mean of the 
correlated signal from adjacent pixels, and subsequently 
removed. The MAMBO pixels are separated in the ar- 
ray by two beams and weak point sources (<150 mJy) 
will therefore not be detected in adjacent pixels. Us- 
ing the measiiremcnts of the atmospheric opacity from 
the skydips, and the elevation at which objects wc;re ob- 
served, the extinction due to atmospheric water vapour is 
then corrected. Gain calibration converts counts to flux 
density, and finally, for each science target the weighted 
mean of all the scans is taken. These are the flux den- 
sities quoted in Table 1: a relatively uniform root mean 
square (rms) noise ~0.55 mJy beam"' was achieved for 
most of the sources, with only a few sources having an 
rms ^0. 7-0.9 mJy beam"'. 

In September 2008 IRAM issued a warning to 
MAMBO observers: several technical problems had 
ocurred during the Win07 semester, which included ther- 
mal and electric malfunctions in the preamplifler, height- 

11 see http://iram.fr/IRAMFR/ARN/dec05/node9.html for 
more details. 



cned sensitivity of the bolometer to microphotonics and 
spikes in the detected signal. These problems were ex- 
pected to affect a fraction of the scans, and the warning 
suggested that any <10 mJy detections obtained dur- 
ing this semester should be checked. Four of sources 
from this sample were detected during Win07: AMS07, 
AMS08 AMS12 and AMS17. All four were observed 
again during the Win08 period, together with extra scans 
of AMSIO, AMSII and AMS20. The data for Win07 
and Win08 were reduced independently, with the results 
shown in Table 2. 

The Win07 flux densities are systematically higher, so 
only sources detected in the Sum08/Win08 period are 
considered safe detections. However, the source AMS17 
was independently detected by Sajina et al. (2008, their 
source 22558), so their flux density is used in Table 1. 
The source AMS07 is not considered a safe detection. 
The potential worry was detecting sources spuriously so 
for the sources that are clearly not detected the combined 
data can be used to yield deeper limits. Table 2 marks 
in bold the flux densities that are used in Table 1. 

2.2. SHARC observations 

Two objects were observed in continuum at 350 /im, 
AMS13 and AMS16, on June 16th and 17th 2005, using 
the Submillimeter High Angular Resolution Camera II 
(SHARC-II) at the 10.4 m telescope of the Cahcch Sub- 
millimeter Observatory (CSO) at Mauna Kea, Hawaii 
(Dowell et al. 2003). SHARC-II is a background-limited 
camera utilizing a "CCD-style" bolometer array with 12 
X 32 pixels. At 350 /im, the beam size is 8.5 arcsec, 
with a 2.59 x 0.97 arcmin^ field of view. The observa- 
tions were taken during very dry weather: AMS16 was 
observed during 2 hours with an opacity at 225 GHz of 
7'225 GHz =0.05, and AMS13 was observed for 1 hour with 
T225 GHz =0.06. The Dish Surface Optimization System 
(DSOS, Leong 2005) was used during the observation 
to correct the dish surface figure for imperfections and 
gravitational deformations as the dish moves in elevation. 

The raw data were reduced with the Comprehen- 
sive Reduction Utflity for SHARG-II, CRUSH - version 
1.40a9-2 (Kovacs 2006). We used the sweep mode of 
SHARC-II to observe the sources: in this mode the tele- 
scope moves in a Lissajous pattern that keeps the central 
regions of the maps fully sampled, but causes the edges 
to be much noisier than the central regions. Pointing was 
checked with Uranus and secondary object CRL_16293- 
2422 several times during the observations. The pointing 
accuracy is better than 3 arcsec. Uranus was also used 
as the flux calibrator. We used Starlink's "stats" pack- 
age to measure the flux of the calibrator and estimate 
the noise of the image. The resulting root-mean-square 
noise values for AMS13 and AMS16 are 7.0 mJy beam"' 
and 5.5 mJy beam^', respectively, and neither source is 
detected at the 3cr level. 

2.3. PdBI observations 

Two sources, AMS12 and AMS16 were observed us- 
ing the Plateau de Bure Interferometer (PdBI) to search 
for molecular gas in the host galaxies, by using carbon 
monoxide (CO) as a tracer of the molecular gas. 

AMS16 was the first object observed with MAMBO, 
in 2005. At the redshift of the source, z = 4.169, the CO 
(4-3) rotational transition falls in the 3 mm atmospheric 
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window [where (4-3) stands for the transition from the 
J = 4 to the J = 3 angular quantum numbers] . 

The observations were carried out on 28 August and 01 
September 2007, with 5 antennae in the most compact 
configuration (D). The source 3C 345 was used as am- 
pHtude and bandpass cahbrator, MWC 349 as amphtude 
cahbrator and the phase cahbrators were 1823+568 (28 
Aug.) and 1637+574 (01 Sep.). The correlator was con- 
figured with both polarizations observing at a central fre- 
quency of 89.2002 GHz, with a bandwidth of 1 GHz. For 
the CO (4-3) transition at 2; = 4.169, this corresponds to 
±1680 km s^^, and should suffice to avoid missing the 
line given the accuracy of the redshift. 

The night of 28 August had moderately good condi- 
tions, and 25% of the data were flagged automatically^^. 
The night of 01 September had excellent conditions (with 
only 1% of the data flagged) although for one antenna, 
1 spectral unit out of the 8 (L07 in horizontal polariza- 
tion for antenna 6) had to be flagged for the entire night. 
The total usable time from both nights amounted to 12.9 
hours on-source. 

The data were reduced using the GILD AS software^'^. 
After amplitude, phase and bandpass calibration, data 
cubes of the phase calibrators and of AMS16 were cre- 
ated. The phase calibrators came out perfectly centred 
at the phase center, suggesting excellent phase calibra- 
tion, but AMS16 was not detected either in continuum 
or in line. A spectrum was extracted at the position 
of the source, with the shape of the beam. The root- 
mean square noise achieved in the spectrum was 0.65 mJy 
beam~^ per 10 MHz bin. Assuming a box-car CO line 
shape with full width to zero intensity of 400 km s"-'^, 
this corresponds to 74 mJy beam~^ km s^^. The contin- 
uum sensitivity achieved over the entire bandwidth was 
0.065 mJy beam~^ rms. 

AMS12 was also observed using the PdBI, to look for 
the CO (3-2) transition. The observations, carried out 
during the 30 April and 13 May 2009, were centred on 
91.796 GHz using a 1 GHz bandwidth. Depending on 
the day, the calibrators 3C 345, MWC 349, 2145+067 
and 0923+392 were used as flux or bandpass calibrators, 
while the source 1637+574 was used as a phase calibrator 
on both dates. On the 30 April, only 5 antennae were 
available and in addition approximately 50% of the data 
were flagged. For both dates the water vapour radiome- 
ter for antenna 2 did not work, so no atmospheric phase 
corrections were applied to the data from this antenna. 
In total, 7.3 hours of usable data were obtained and an 
rms noise of 0.70 mJy beam~^ per 10 MHz channel was 
reached. 

The CO (3-2) transition was clearly detected in 
AMS12, with an integrated flux of 630 mJy beam~^ 
km s^^. The continuum rms is 0.1 mJy beam^^, and 
no continuum is detected at 3 mm, so that we obtain a 
3(7 limit of 5*3 mm <0.3 mJy. In the case of AMS16, the 
lack of continuum sets an Sa upper limit on the 89 GHz 
(3 mm) flux density of S3 mm <0.2 mJy. 

Data can be flagged as having bad quality for many reasons. 
Examples for reasons to automatically flag the data during the 
observations are time errors, tracking problems or antenna shad- 
owing. During the calibration and reduction of the data, visibilities 
with high phase noise or high pointing correction, tracking errors 
as well as amplitude loss can also be flagged. 
13 http: //www.iram.fr/IRAMFR/GILDAS 



3. RESULTS FROM CONTINUUM OBSERVATIONS 

The millimetre continuum observations provide us with 

a wealth of information about the far-infrared luminosi- 
ties and cool-dust masses of our sources. They can also 
be used to test whether there is enough dust along the 
host galaxy to cause the obscuration in the quasars. In 
the following section, we derive these physical properties 
from our observations. 

The data uniformly available for this sample includes 
the following wavelengths: 3.6, 4.5, 5.8, 8.0, 24, 70 and 
160 yttm, as well as 1.2 mm. For a small number of sources 
data is also available at 350 or 850 /^m. With such ex- 
cellent wavelength coverage, one might expect the far- 
infrared luminosities and cool-dust masses to be very 
accurately constrained. However, the wavelengths cor- 
responding to the rest- frame mid-infrared arc dominated 
by hot dust from the warm torus around the central en- 
gine, and cannot be used to constrain the cool-dust com- 
ponent as we explain below (this was originally by selec- 
tion, but has additionally been confirmed by mid-infrared 
spectroscopy, see Martinez-Sansigre et al. 2006b, 2008). 

The far-infrared luminosity is dominated by cool dust, 
expected to be characterised by temperatures ~35-50 K. 
This means that emission is expected to peak around 65- 
90 fim and drop exponentially at shorter wavelengths. 
The emission from the warm torus is caused by dust at 
a range of temperatures, typically ~100-2000 K, whose 
emission will peak at shorter wavelengths than that of 
the cool dust. If the AGN is powerful enough, as is the 
case in the sources making up this sample, the dust emis- 
sion from the torus will dominate at mid-infrared wave- 
lengths. Hence, the mid-infrared emission is not useful 
in constraining the cool dust mass and luminosity. 

To illustrate this. Figure 1 shows a model rest-frame 
SED consisiting of cool dust modelled by a gray body 
with temperature T = 40 K and emissivity (3—1.5 to- 
gether with a model AGN torus from the radiative trans- 
fer models of Nenkova et al. (2008a,b). The Spitzer MIPS 
bands at 24, 70 and 160 /xm, as well as the MAMBO 
band at 1.2 mm have been marked in gray, assuming 
this source to be at z = 2. The emission seen at 1.2 mm 
will be completely dominated by the cool dust, while the 
emission at 24 /xm will be completely dominated by the 
warm dust of the torus. The 70 /zm band is expected 
to be dominated by the torus, with some contribution 
from the cooler dust, while the situation for the emission 
observed at 160 /zm is expected to be reverse: the cool 
dust dominates with some contribution from the torus. 

This is only an illustration, and the relative impor- 
tance of each component will vary depending on the ac- 
tual temperature of the cool dust, the actual SED of the 
torus and their relative luminosities. However, Figure 1 
does show that when a powerful AGN is present, the 
data at 24 and 70 /zm and the mid-infrared spectra will 
not provide very powerful constraints on the cool-dust 
emission. It is the data at 1.2 mm that will provide the 
best measurement of this emission. The data at 160 /im 
should provide useful constraints but for most sources in 
our sample it is too shallow to provide detections. Data 
at 350 or 850 /xm is extremely useful but is only available 
for 3 sources. 

Hence, we will determine the far-infrared luminosities 
based on the flux densities at 1.2 mm alone for most of 
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the sources in this sample. For sources where the data 
at 160, 350 or 850 fim provide a useful constrain, we will 
make use of this data too. However, the 24 /xm data will 
not be used to constrain the cool dust emission (and in 
most cases, neither will the 70 /um data). 

3.1. Flux densities at 1.2 mm 

The MAMBO observations yield detections at the >3a 
level for 5 out of 21 sources observed (24%). The results 
are summarised in Table 1 and Figure 2 (where sources 
with no spectroscopic redshift were placed at z = 1.98, 
the median for sources with a spectroscopic redshift). 
The mean noise achieved is 0.56 mJy bearn"^, the me- 
dian is 0.54 mJy beam~^. Hence, for the non detections, 
typical 3(7 limits ^^1.65 mJy are obtained. 

With 76% of the sources undetected, the median flux 
density does not provide a useful characteristic number 
for the sample. A more meaningful number is the mean 
flux density at 1.2 mm for the entire sample. The flux 
densities of all sources were therefore combined using a 
weighted mean: 




where Wi = l/af. This yielded a statistical detection, 
with (5i.2 mm) = 0.96 ± 0.11 mJy. A potential worry is 
that a few bright sources might be causing the statistical 
detection, while the rest of the sources arc intrinsically 
faint. Stacking only the non-detections still yields a de- 
tection at the 3.9cr level: (51.2 mm) = 0.51 ± 0.13 mJy. 
Stacking all the detections leads to a mean flirx density of 
{S1.2 mm) = 2.30 ± 0.23 mJy. Figure 3 shows the distri- 
bution of flux densities (from Table 1), which is centred 
on a positive value and skewed towards the higher flux 
densities. This flgure supports the statement that the 
typical flux density of an object from this sample is ~0.5- 
1.0 mJy. If the actual flux densities of the sources were 
so low that the distribution were dominated by the noise. 
Figure 3 would be centred on 0, with a symmetric charac- 
teristic width ~0.55 mJy, and only a slight skew towards 
the high fluxes (the detections with Sv<^1.65 mJy). 

We now consider whether the detection at 1.2 mm is 
related to the presence or lack of narrow lines in the 
optical spectra. Of the 5 detections at 1.2 mm, 4 show 
narrow lines (80%) and only 1 (20%) has a blank optical 
spectrum (AMS19). Alternatively, of the 9 objects with 
narrow lines, 4 have detections at 1.2 mm (44%), while 
of the 12 with no narrow lines, only 1 is detected (8%). 
Due to the small numbers, no trend is significant, but at 
first sight the data hints that sources with narrow lines 
are more likely to be detected at 1.2 mm. 

Looking at the stacked flux densities for narrow-line 
and for blank objects (no narrow lines) suggests the dust 
content to be similar: the weighted mean for narrow- 
line sources is 1.09±0.16 mJy, while for blank sources 
it is 0.81±0.16 mJy, so both are statistically detected. 
Although at first sight the difference appears significant, 
the narrow-line objects include AMS12 with a fiux of 
3.7 mJy which is affecting the mean: the mean of the 



other 8 narrow line objects excluding AMS12 is 0.76 mJy, 
very similar to that of the objects with no narrow lines. 
Hence, the mean mm flux densities for both narrow line 
and blank objects appear similar. 

3.2. Far-infrared luminosities 

Given a flux density, the infrared luminosity of a gray 
body is given by: 

^FIR = 

15 Andj f kB y r(4 + ,j)C(4 + J)asBr^+'5i.2nu^ 

where hp is Planck's constant, fee is Boltzmann's con- 
stant, (TsB is the Stefan-Boltzmann constant, and c^l is 
the luminosity distance at redshift z. The terms v and 
frf refer to observed and rest-frame frequency respec- 
tively, Urf = (1 + z)u. The term B^,^^ (T, v^f) is the Planck 
radiation function: 

Bu^, (T, z/rf' ) = — ^ \ ^ ' 

I e"^ - 1 J 

Figure 4 and Table 3 show the inferred values of Lpm 
assuming two different gray bodies, one with T = 47 K 
and P — (representative of unobscured quasars at 
high redshift, see Beelen et al. 2006), the other with 
T = 35 K and /3 = 1.5 (representative of submillimetre- 
selected galaxies, SMGs, see Kovacs et al. 2006). The 
value of (5*1 2 mm), also assuming z = 1.98, corresponds 
to (Lfir) =6.3±0. 7x10^2 (r = 47 K and /? = 1.6) 
or 1.6±0. 2x10^2 Lq (for T = 35 K and /3 = 1.5). 
Hence, the obscured quasars are typically ultra-luminous 
in the far-infrared, with a characteristic luminosity of 
~ 4 X 10""^^ L0 (the mean of both values). 

Despite the range in spectroscopic redshifts (1.6 < 
z <4.2), the relatively flat selection function at 1.2 mm 
means that the detection of sources is primarily de- 
termined by their far-infrared luminosities, and is only 
weakly dependent on redshift (see Figure 4). In fact, for 
a given flux density, the corresponding far-infrared lumi- 
nosity is actually smaller at high redshift. Hence intrin- 
sically less luminous objects are easier to detect at higher 
redshift. This probably explains why in our sample, the 
higher redshift sources (z ^2) have a slightly higher de- 
tection rate than the lower redshift sources [z ^2). 

The luminosity arising in the far-infrared is a consid- 
erable fraction of the AGN bolometric luminosity (Lbob 
between infrared and X-ray energies, so between 10^ and 
10""^^ Hz) which is typically 2 x 10^'^ for our sample 
(see Section 4). Hence, Lfir is typically ~0.1-0.3Lboi 
depending on the assumed gray body. This is slightly 
higher than what is found in low-redshift unobscured 
quasars, where Z/piR^O.liboi (e.g. Sanders et al. 1989; 
Elvis et al. 1994; Schweitzer et al. 2006, and see also 
Section 7). For SMGs, LpiR^Lboi since the far-infrared 
luminosity consitutes the bulk of the bolometric emission 
(e.g. Egami et al. 2004). 

In addition to the data at 1.2 mm, we have access to 
data at 70 and 160 /im for all 21 sources (from Frayer 
et al. 2006), data at 350 for 2 sources (AMS13 and 
AMS16, see Section 2.2) and data at 850 /xm for one 
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source (AMSll, data from Frayer et al. 2004). In most 
of the cases, the two fiducial gray bodies provide an ac- 
ceptable fit (see Figure 5). For the sources where neither 
the T = 47 K and /? = 1.6 or the T = 35 K and /? = 1.5 
gray bodies are appropriate, we use the data to find bet- 
ter fits. We only fit sources that have been detected at 
1.2 mm. 

For AMS08, Figure 5 suggests that the 70 fiui data 
are dominated by the cool dust component and not by 
the warm torus, so we use the data at 70 and 160 /xm, 
as well as 1.2 mm, and allow (3, T and Lfir to vary. 
For AMS12 and AMS19, the 70 fim data do not provide 
tight constraints, so /3 is kept fixed at a value of 1.5, and 
only T and £fir arc allowed to vary. The results are 
summarised in Table 4 and overplotted in Figure 5. 

3.3. Estimated star-formation rates 

Under the assumption that the far-infrared luminosity 
arises solely from dust heated by young stars, and not by 
the AGN itself, star-formation rates (SFRs) can be esti- 
mated from LpiR- Since the dust is primarily heated by 
young stars, the conversion to total SFR involves assum- 
ing an initial mass function (IMF) . We use the conversion 
from Kennicutt (1998), 

where Xir, refers to the total infrared luminosity. This 
conversion assumes a Salpeter (1955) IMF. Given that at 
mid-infrared wavelengths our sources are dominated by 
emission from warm dust (T >100 K) heated by the AGN 
and not by stars (see Martfnez-Sansigre et al. 2008), 
the mid-infrared luminosities due to star-formation can- 
not be easily derived. We therefore approximate the 
total infrared luminosity to the far-infrared luminosity, 
^IR~^FIR,. -^FiR is determined by integrating analyti- 
cally the fiducial gray-bodies at all wavelengths, so there 
is no wavelength at which the integration of the gray 
body is cut off. However, any other components which 
would contribute to Lir and are not part of the gray 
body (e.g. PAHs or mid-infrared continuum) are ne- 
glected. Table 3 summarises the results. 

For the sources detected at 1.2 mm, the inferred values 
of XpiR and SFRs are in the range 300-3000 M© yr"!, 
comparable to those of SMGs (e.g. Barger et al. 1998; 
Hughes et al. 1998; Chapman et al. 2004; Greve et al. 
2005). The SFR corresponding to (Lfir) - 4 x lO^^ L© 
suggests that the characteristic SFR for the sample is 
~700 M0 yr-i. 

Martmez-Sansigre et al. (2008) detected polycyclic 
aromatic hydrocarbon (PAH) emission in three of the 
sources in this sample (see Table 1), and found the lumi- 
nosity of the detected 7.7 fim PAHs to be comparable to 
the typical luminosities of PAHs found in the sample of 
SMGs studied by Valiante et al. (2007). The far-infrared 
luminosities inferred here are slightly lower but still com- 
parable to those of SMGs. 

If a significant fraction of Lfir is actually due to AGN- 
heated dust, these values for the SFRs will be overesti- 
mated. This could happen if dust at i^kpc scales is being 
heated by the Lboi^ lO^'' Lq quasar (see e.g. Barvainis 
1987; Sanders et al. 1989; Siebenmorgen et al. 2004, and 
Section 3.5). 



In addition to the infrared data, radio observations are 
also available, which could in principle yield further con- 
straints on the SFRs. The multi-frequency and high- 
resolution radio data, however, suggest the radio emis- 
sion is completely AGN-dominated, and therefore does 
not provide information on the SFRs (see Martfnez- 
Sansigre et al. 2006a; Klockner et al. 2009). 

It is not possible to discriminate whether the far- 
infrared luminosity arises from AGN- or star-formation- 
heated dust: PAHs were only detected in 3 objects in our 
sample yet the upper limits obtained for the other sources 
correspond to very high PAH luminosities (and hence are 
not very stringent limits, see Martfnez-Sansigre et al. 
2008). The SFRs inferred for our sample are extremely 
high but physically-possible. 

3.4. Cool-dust masses 

The mass of cool dust responsible for the far-infrared 
emission can be estimated as: 

71^ _ d1Si,2 mm 

{l-\-z)Kiiy,{)B,^,iT,u,{)' ^' 

where k is the mass absorption coefEcent and the main 
source of imccrtainty, so we will consider two different 
fiducial values to illustrate the estimated values of M^- 
We follow Beclcn et al. (2006) and assume k =0.04 m^ 
kg~^ at 1.2 mm (within the values quoted by Alton 
et al. 2004). The values for individual objects are listed 
in Table 3. Using (51.2 mm) = 0.96 ± 0.11 and again 
assuming z = 1.98 as the typical redshift, we estimate 
the characteristic dust mass to be {M^)— 2.9 x 10* M© 
and 5.1x10* M© (for T = 47 K, /3 = 1.6 and T = 35 K, 
/3 = 1.5 respectively). If instead we had followed Dunne 
et al. (2003) in assuming a value of k =2.64 m^ kg"-'^ 
at 125 /im, the resulting estimates of (Md) would have 
been = 1.6 x 10* Mq and 2.3x10* Mo(for T = 47 K, 
/3 = 1.6 and for T = 35 K, /3 = 1.5), which makes little 
difference to our arguments. In any case, the values of 
Md quoted should be taken as indicative only, and we 
take 3 x 10* Mq as the fiducial cool-dust mass for this 
sample (the mean of the four values quoted above). 

3.5. Characteristic scale of cool dust 

We now estimate the characteristic distance of the dust 
under two assumptions: that it is heated by young stars, 
and that it is heated by the AGN. 

We note that if the dust is heated by stars, then the 
inferred SFRs are similar to those of SMGs. Interfer- 
ometric observations of SMGs suggest a characteristic 
scale of ~2 kpc (Greve et al. 2005; Tacconi et al. 2006; 
Younger et al. 2008), so if the cool dust around our ob- 
scured quasars is heated by stars, the scale of the dust is 
expected to be similar to this value. 

To estimate the scale of AGN- heated dust, wc cqiiate 
the grain rate absorption of energy from ultraviolet pho- 
tons to the radiated energy (see e.g. Barvainis 1987, for 
a derivation): 

(6) 
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where L^v is the uhraviolet luminosity of the central 
source, and Tuv is the opacity at ultraviolet wavelengths. 
The symbols a and p stand for the grain size and mate- 

rial density, and Q{a, v) = / j is the emissivity. 

Following Alton et al. (2004) we assume a =0.1 /xm and 
p =3000 kg m~^ as fiducial values, which yields: 



(kpc) 



16.3 



T 
40 K 



5 X 1012 L, 



(7) 



The typical extinction-corrected bolometric luminosity 
of our sample of quasars is 2 x IQ^'^ Lq (see Section 4). 
Assuming 

-^uv ^ 0.25 -^boi (^'S' Elvis et al. 1994) the 
typical Z/uv is about 5x10^^ Lq. If we assume the ul- 
traviolet photons travel unhindered up to the character- 
istic radius of the dust, so that = 0, and assume 
T ~ 47 K and /? = 1.6, then the dust responsible for 
the 1.2 mm emission is at a distance from the AGN of 
~10 kpc (27 kpc for T = 35 K and /3 = 1.5). 

Given that t^v is likely to be > at kpc scales, since the 
cool dust itself can absorb the uv photons from the AGN 
and shield the dust on larger scales (see Section 3.6), 
these scales should be taken as upper limits. Hence the 
cool dust is expected to be distributed on scales ;^10 kpc. 

3.6. Host obscuration? 

The characteristic cool-dust mass derived from our ob- 
servations can also be used to estimate how much obscu- 
ration the host galaxy can cause. We are considering the 
cool dust only {T Si50 K), which emits in the far- infrared. 
We are interested in the amount of extinction it can cause 
at optical and mid-infrared wavelengths, where the emis- 
sion from the cool dust is negligible. Hence, a simple 
screen approximation neglecting the emission from the 
obscuring dust is appropriate, so that a full radiative- 
transfer calculation can be avoided. 

In order to do this, we use the dust extinction models 
of Pei (1992) for the Milky Way (MW) and Small Magel- 
lanic Cloud (SMC). These laws, however, correspond to 
/? = 2.0 for A <l200 /im, so extrapolating directly from n 
at 1.2 mm to the visual band (0.55 /xm) using these ex- 
tinction models would not be consistent with the values 
of/3 = 1.5 — 1.6we have assumed. 

We therefore estimate tlie opacity at 200 /im, given a 
typical dust mass of 3x10® Mq and assuming /3 = 1.6. 
Assuming a spherically-symmetric smooth distribution 
of dust with a characteristic radius r and volume V = 
the column density and hence the oppacity can be 
estimated: 



T200 /nm = r-^K200 /im- 



(8) 



Then, using an extinction law, the corresponding ex- 
tinction in the visual band. Ay, can be derived (in 
units of magnitudes). Assuming a characteristic radius 
of 10 kpc, K =0.04 m^ kg~^ at 1.2 mm, we estimate 
^200 urn = 1-2 x 10~^, which requires Ay='l for the MW 
model and 1.5 for SMC. If instead we assume k =2.64 m^ 
kg~i at 125 /im and the same radius, the opacity in- 
creases to T200 /im = 2 X lO^'^ and requires Ay=7.5 (MW) 
or 2.5 (SMC). 



Assuming instead a radius of 2 kpc, characteristic of 
the gas and dust in SMGs (Greve et al. 2005; Tacconi 
et al. 2006; Younger et al. 2008), the inferred value of 
T200 /im is 0.03 — 0.05, corresponding to values of Av=35- 
65 (SMC dust and depending on k) or even > 100 (MW). 
A radius of 20 kpc leads instead to low values of Ay 
~0.3-1.5 (using both MW and SMC extinction curves). 

Given that the condition for a quasar to be considered 
obscured is Ay>5 (see Simpson et al. 1999), and com- 
paring to the estimates of Ay from Section 4 shown in 
Table 3, the main conclusion from this estimate is that 
there is probably enough cool dust in the host galaxies of 
many of these sources to obscure the quasars, regardless 
of the orientation of the torus with our line of sight. 

There are, however, examples of high-redshift unob- 
scured quasars with similar (sub)millimeter fluxes and 
hence similarly large dust masses (e.g. Omont et al. 2003; 
Beelen et al. 2006): the presence of large dust masses 
does not guarantee heavy obscuration. This is also the 
case for our sample (see Section 3.1). A likely explana- 
tion for this is that the dust in the host galaxy has a 
very irregular or clumpy distribution, so that only some 
lines of sight will be obscured. Rather than orientation- 
dependent obscuration, as is the case for orientation 
caused by the torus, the obscuration (or lack of) by dust 
in the host galaxy is likely to be dependent on luck, al- 
though samples of obscured quasars are expected to be 
biased towards "bad luck" in the form of blocked lines of 
sight. 

4. MID-INFRARED SED FITTING 

In this section, we make use of broad-band data to 
characterise the SEDs of our sources. We use the IRAC 
and MIPS-24 /xm data to obtain an estimate of Lboi 
and from broad-band SED fitting. The approxi- 
mation we use is that an obscured quasar SED can be 
parametrised as an unobscured quasar SED, with bolo- 
metric luminosity Lboi, with a foreground screen of dust 
causing an amount Ay of extinction. 

The values for Xboi and Ay obtained in this way are 
only a parametrisation of the obscured SED. The Ay is 
only the apparent extinction, assuming the unobscured 
SED to be the underlying SED. This is an accurate de- 
scription at optical and ultraviolet wavelengths: the in- 
trinsic illuminating source is the accretion disk around 
the SMBH, which emits thermally with a range of tem- 
peratures ^ 10^ — 10^ K, and the dust will simply absorb 
this and re-emit it at longer wavelenghts. 

The mid-infrared emission from unobscured quasars, 
however, is already reprocessed light from warm dust 
(T ~100-2000 K). The dust responsible for this emis- 
sion must be directly illuminated by the accretion disk, 
and must be relatively close to the source of ultraviolet 
photons: it is likely to originate in the inner region of 
the torus. For example, following Section 3.5 and using 
Equation 7, T ^300 K corresponds to ~50 pc. For an 
unobscured quasar, the mid-infrared emission will travel 
almost unhindered until it reaches the observer, while in 
the case of an obscured quasar, it still has to go through 
a considerable amount of dust. Most of this dust will be 
further away (and thus colder) and might not be directly 
illuminated by the central engine anyway. Thus, the dust 
obscuring the mid-infrared emission is once again un- 
likely to emit very significantly at the same wavelengths. 
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The approximation of the obscured SED as an unob- 
scured SED with a foreground screen of dust is still de- 
fendable at mid- infrared wavelengths, although it will be 
a poorer approximation than in the optical regime. In 
any case, it is a simplification of the situation. 

This approximation breaks down completely in the far- 
infrared, so we do not attempt to model the far-infrared 
SEDs in this way. Instead, the far-infrared component 
of the SEDs was modelled as a gray body, as described 
in Section 3.2. 

To model the SEDs in this way, we use the 24-/zm data 
from the catalogue of Fadda ct al. (2006), and the IRAC 
data from Lacy et al. (2005) as well as deeper observa- 
tions around 4 objects from our program PID 20705 (PI 
M. Lacy, the objects are AMS05, AMS12, AMS16 and 
AMS17). We use the Elvis et al. (1994) median quasar 
SED. The quasar SED is then obscured by a screen of 
MW-type dust, using the models of Pei (1992). It is this 
dust law that will determine the depth of the silicate fea- 
ture (around 9.7 iJ.ni) in the models. An elliptical galaxy 
SED (from Coleman et al. 1980) is also used, to repre- 
sent the old stellar population, and this stellar light is 
not subject to any obscuration. The quasar bolometric 
luminosity. Ay, and the luminosity of the galaxy are all 
allowed to vary, and the best fit is kept^'*. 

For sources with many limits (namely AMS08, AMSIO, 
AMS15, AMS18), the marginalised posterior distribution 
function (PDF) was flat for a range of values of Ay (i.e. 
the same ^^.s obtained for all values of Ay above 
a certain value). Given that low-luminosity objects are 
always more common than high-luminosity objects, and 
given that the inferred value of Lboi will correlate with 
Ay^^, we chose the lowest value of Ay (to the closest 
multiple of 5) and associated Lboi within the flat region 
of the posterior PDF. The results are summarised in Ta- 
ble 3 and all SEDs are presented in Figure 5. The best-fit 
SEDs are shown (blue solid line) , as well as the individual 
components (quasar and elliptical galaxy, dashed lines). 
The broad-band data at 3.6, 4.5, 5.8, 8.0, 24, 70 and 
160 nm as well as the IRS spectra are superimposed. 
On the far-infrared end, the 1.2 mm point, and the two 
gray bodies are also shown. For the sources AMS13 and 
AMS16, the 350 /im data from SHARC are also shown, 
and for AMSll a data point at 850 //m from Frayer et al. 
(2004) is also included. For source AMS05, additional 
data at 1.2, 1.7 and 2.2 /xm are also shown (from Smith 
et al. 2009). 

The IRS spectra were not used in the fit, but used 
to provide instead an independent quality check on the 
fits. In most sources the best fit to the broad-band SED 
is in reasonable agreement with the IRS spectrum. In a 
few cases the IRS spectra show deeper silicate absorption 

^■^ This is a essentially the same routine as used in Martmez- 
Sansigre et al. (2007), except the redshift is not varied and their 
"blue" component is not included: we use the spectroscopic red- 
shifts, or assume z = 1.98 when no redshift is available, while 
the blue component is irrelevant in the wavelength range discussed 
here (3.6-24 /im). In addition, we are not selecting models here, 
only fitting parameters, so for each object we simply keep the val- 
ues with the lowest x^- Finally, due to the narrower wavelength 
range here, we allow the errors in photometry to be < 10% (c.f. 
Martmez-Sansigre et al. 2007). 

In order to match the observed flux density of a given source 
at a given redshift, increasing the Ay will cause the inferred value 
of Lboi to increase. 



features than those resulting from the best-fit Ay (e.g. 
AMS06, AMS19), in other cases the agreement is very 
good (e.g. AMS05, AMS13). The overall agreement be- 
tween the IRS spectra and the best-fit broad-band SEDs 
is reasonably good. 

Two sources have very low values of Ay from the 
broad-band SED fitting: AMS02 (Av=0.5) and AMS20 
(^v=0.0). This would suggest both objects are unob- 
scured quasars, and should not make it into the sam- 
ple. However, both objects showed blank optical spec- 
tra, despite long integrations (see Martfnez-Sansigre 
et al. 2005, 2006b). The IRS spectrum of AMS02 shows 
the silicate feature in absorption while the spectrum of 
AMS20 shows a hint of silicate absorption (although in 
Martfnez-Sansigre et al. 2008, it was not considered se- 
cure due to the low signal-to-noise ratio). Both sources 
have 3.6 /im flux densities at the edge of the selection cri- 
teria, and AMS02 also has a 24 /im density at the edge 
of the criteria (see Table 1 and Figure 1 of Martfnez- 
Sansigre et al. 2005). AMS02 is certainly a heavily 
obscured source, the flat SED might be due to a rela- 
tively bright host galaxy compared to the AGN emission. 
AMS20 is less clear, the flat SED could also be due to 
the host galaxy contribution, and in any case the blank 
optical spectrum shows it is certainly not an unobscured 
quasar. It might be a reddened quasar {0<Ay^5), or an 
obscured quasar with a pecifliar mid-infrared SED. 

For our sample, we find no significant correlations be- 
tween Xfir and iboi or Lfir and Ay. 

5. A MEAN OBSCURED QUASAR SED 

We now derive a typical dust SED for our sample of 
obscured quasars and compare it to clumpy models of 
AGN tori. 

As mentioned earlier, due to the large number of up- 
per limits the median flux density at 1.2 mm is not a 

very useful quantity (since it yields an upper limit), yet 
we find the weighted mean yields a statistical detection. 
To be consistent, we will therefore consider mean quan- 
tities for the mid-infrared too. We note that within our 
sample there is no correlation between iboi and Xfir, so 
wc treat the (rest-frame) near-infrared, mid-infrared and 
far-infrared components independently. The mean SED 
of our spectroscopically-confirmed objects is shown as a 
solid red line in Figure 6. 

We do warn about several caveats: The far-infrared, 
mid-infrared and near-infrared components have been 
determined in different ways. The far-infrared compo- 
nent is an analytic gray body to the mean flux density 
at 1.2 mm, with an empirical quasar SED (from the 
mean Z/boi)- It is still ill constrained, since only some 
objects have been detected at 70 or 160 /im. The mid- 
infrared component results from stacking IRS spectra, 
and is probably the best determined part, although it 
uses only sources with spectroscopic redshifts and IRS 
spectra. 

For the entire SED, only sources with spectroscopic 
redshifts are used. For sources with no redshift from op- 
tical spectroscopy, redshifts can only be determined if the 
silicate feature is deep enough to be noticeable in noisy 
spectra. The mean mid-infrared SED is therefore par- 
tially biased towards sources with deep absorption fea- 
tures at 9.7 /xm. 

The near-infrared component uses the data from all 
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sources, but because these sources have a variety of red- 
shifts, of values of Ay, and of host galaxy luminosities, 
it is really a mixture of heavily-extinct AGN light and 
stellar light. 

In addition, the sources in this sample probably in- 
clude examples of both torus- and host-obscured quasars, 
meaning that the mean SED will represent a mixture of 
such objects. Despite all these caveats, we consider the 
mean SED to be a useful fiducial representation of ob- 
scured quasars. 

5.1. Far-infrared component 

We parametrise the far-infrared component as being 
described by the two fiducial gray bodies. The far- 
infrared luminosity is determined by the weighted mean 
flux density, {S1.2 mm) = 0.96 ± 0.11 mJy at z = 1.98. 
Assuming the gray body with T =47 K and /? =1.6 re- 
sults in (Lfir) = 6.3 X 10^^ Lq, while assuming T =35 K 
and P =1.5, we obtain (Lfir) = 1-6 x 10^^ Lq. 

5.2. Mid-infrared component 

The mid-infrared component was obtained by shifting 

the IRS spectra to rest-frame wavelengths, interpolat- 
ing to a common wavelength grid and taking the mean 
luminosity density of the 14 sources with spectroscopic 
redshifts and mid-infrared spectra. Different wavelength 
bins have a different total number of sources contribut- 
ing, and it was decided to cut below 4.5 /um and above 
12.5 /txm, where 3 or less sources were contributing. Fig- 
ure 6 shows the mean IRS spectrum smoothed by a box- 
car average of 0.7 /zm. 

The depth of the silicate absorption feature at 9.7 fim, 
T9,r is usually measured by interpolating the continuum 
at 9.7 ^m using the continuum bluewards and redwards. 
The difference in flux densities between the interpolated 
continuum (/ic) and the observed feature (/of) then 
yields rg.y: 

r„.l„(|^) (9) 

This is difficult to measure from the mean IRS spec- 
trum, since there are no data points redwards of the ab- 
sorption feature, so no real interpolation can be made. 
However, comparing the flux densities of the 9.7 fxm fea- 
ture and the 8 /zm continuum in Figure 6, the depth of 
the feature is likely to be in the range 2.5 ^ T9.7 ^ 3.0. 

Figure 6 also shows that this is deeper than what is 
expected from the unobscured SED with the screen of 
dust (using the mean Ay. dashed line), which otherwise 
fits the region of 4.5-8.5 /xm quite well. We note that the 
mean Ay from all sources with spectroscopic redshifts 
(j4v=32 is very similar to the mean of the 14 sources used 
for the mid-infrared component. Ay =34). This small 
difference leads to a practically indistinguishable SED in 
the mid-infrared and cannot explain the discrepancy in 
the depth of the silicate feature. 

At this stage, the discrepancy with the silicate feature 
is only due to the dust extinction model used: an ex- 
tinction law with deeper silicate absorption (such as the 
small magellanic cloud dust model of Pei et al. 1992) can 
plausibly reproduce the observed feature. 

However, a typical silicate feature as deep as ^ 2.5 
suggests some degree of geometrical obscuration caused 



by cold dust (see Levenson et al. 2007, for a discussion). 
This lends further support for some degree of host ob- 
scuration. 

Finally, the mean spectrum from the IRS shows an ab- 
sorption feature at ~7 /xm probably due to hydrogenated 
amorphous carbons (see e.g. Spoon et al. 2002). 

5.3. Near-infrared component 

For the near-infrared component, for each object the 
flux densities of the IRAC bands were converted to lu- 
minosities at the rest-frame wavelength. They were then 
binned in a grid ranging between 0.5 and 4.5 /im. Look- 
ing at the SEDs in Figure 5 it seems that the upper limits 
are generally close to the detections in adjacent bands, so 
we treat limits like detections. In the bins with contribu- 
tions from several sources, the mean luminosity density 
was computed. However, some of the bins have a contri- 
bution from only one object and some from none (these 
last bins are not shown in the plots). The solid red line in 
Figure 6 shows the mean photometry at rest-frame near- 
infrared, smoothed by a box-car average of 0.6 /xm. The 
line is irregular, because of the vastly different individual 
SEDs seen in Figure 5: in some objects, all IRAC bands 
are likely to be dominated by the AGN (e.g. AMS19), 
in others by the host galaxy (e.g. AMS16). If instead 
of using limits as data, we do not use them at all, the 
resulting SED is very similar. 

This part of the SED should clearly be treated with 
caution. However, an estimate of the luminosity of the 
typical host galaxy can be attempted. Overplotted on 
Figure 6 is a 2; = elliptical galaxy from Coleman 
et al. (1980), normalised to be 12.6 x brighter than an 
L* galaxy in the local K-band luminosity function (Cole 
et al. 2001). Assuming ~2 magnitudes of passive evo- 
lution between z ^ 2 and z = 0, it is the equivalent to 
a present day 2L* galaxy (this is partly by selection, as 
is described in detail by Martmez-Sansigre et al. 2005, 
2006b). We can thus see that the luminosity of the near- 
infrared component of the typical obscured quasar SED 
corresponds approximately to the progenitor of a 2L* 
galaxy, although it might also have an important contri- 
bution from AGN light (in which case the galax;y lumi- 
nosity will be overestimated). 

5.4. Comparison to clumpy models of the torus 

It is of particular interest to compare the mean SED 

to clumpy models of the AGN torus. This is shown in 
Figure 7 (left panel), where the mean SED is overplotted 
on the models of Nenkova et al. (2008a, b). The models 
cannot reproduce simultaneously the relatively flat SED 
at wavelengths shorter than 8 jum, the sharp drop in 
the near-infrared and the deep silicate absorption. The 
left panel of Figure 7 also shows that the models with 
deeper silicate features still struggle to match the ob- 
served silicate depth (rg.y ^ 2.5) yet will overpredict the 
far-infrared luminosity. 

The right panel of Figure 7 shows the two flattest 
clumpy SEDs from the right panel, now with additional 
extinction of Ay=8i) added due to foreground dust. 

We see that the clumpy models, with the additional 
extinction and re-emission from the cool dust can repro- 
duce reasonably well the observed SED. Since the dust 
is cold, it re-emits in the far-infrared, meaning it does 
not "fill-up" the silicate absorption feature, unlike dust 
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at ^300 K would (again, see Levenson et al. 2007). To 
reproduce the deep silicate feature, a value of j4v~80 is 
required (see Figure 7). 

The typical dust mass of 3 x 10* Mq can cause this 
extinction of ^v=80 provided that the dust is on scales 
<2. 3-3.0 kpc (from Equation 8 in Section 3.6 and assum- 
ing MW dust and the two values of K200 /im)- This is in 
excellent agreement with the inferred scales of the cool 
dust emission in SMGs (typically ~2 kpc Greve et al. 
2005; Tacconi et al. 2006; Younger et al. 2008). 

Polletta et al. (2008) studied in detail the mid-infrared 
spectra and mid- to near- infrared SEDs of 21 high- 
redshift heavily obscured quasars, selected to be bright 
(>1 mJy) at 24 /im, optically faint (typically R > 24) 
and with mid-infrared spectra characteristic of heavily 
obscured AGN (see their Section 2 for more details). 
They found that 9 out of 21 sources did not fit clumpy 
models well and required an additional foreground screen 
of dust that did not re-emit at mid- or near-infrared 
wavelengths. 

Our findings for the mean SED are consistent with 

their results, and in addition, we are able to estimate 
the Ay in a consistent way from Md, which was derived 
from our observations at 1.2 mm. The resulting dust 
size required for to cause the Ay (~2. 3-3.0 kpc) 
is very close to the sizes estimated for the far-infrared 
emission of SMGs. Again, our results provide evidence 
for obscuration due to dust on kpc scales, dust that has 
lower temperatures and is distributed on larger scales 
than what is expected for the dust of the torus. 

5.5. Com,parison to local ULIRGs 

Figure 8 shows the mid-infrared component of the 
mean SED with the spectra of four local ultra-luminous 

infrared galaxies ovcrplottcd (ULIRGs, with Lm > 10"'^^ 
L0). The spectra for these sources were obtained from 
Armus et al. (2007) and Spoon et al. (2007). These 
ULIRGs represent one pure-AGN, IRAS 12514-1-1027 
(hereafter IRAS 12514), two AGN-starburst composites, 
Mrk 231 and Mrk 273, and one starburst-dominatcd 
source, Arp 220. All spectra have been normalised to 
match our mean SED between 7.0 and 7.8 fim (around 
4x1013 Hz). 

Compared to these sources, our mean SED is relatively 
flat, similarly to Mrk 231 or IRAS 12514, yet the silicate 
absorption feature is significantly deeper, more similar 
to that of Mrk 273 or even Arp 220. Once again, we see 
the peculiar spectral shape of this mean SED, where the 
silicate feature is suprisingly deep given the relative flat- 
ness of the continuum. The most similar nearby ULIRG 
is IRAS 12514, which is known to be an obscured quasar 
(Wilman et al. 2003), although this source has a shal- 
lower feature at 9.7 /xm. 

6. RESULTS FROM OBSERVATIONS OF CO 
6.1. Line properties 

Two sources, AMS12 and AMS16, were observed with 
the IRAM PdBI to search for the CO (3-2) or (4-3) ro- 
tational transitions, tracers of molecular gas. 

The CO (3-2) line in AMS12 is clearly detected, and 
is well described by a gaussian with a peak flux of 
2.14 mJy and a full-width half-maximum (FWHM) of 
275 km s~i (see Figure 9). The integrated line flux of 



this gaussian is Sco^v =630±50 mJy km s^^, where 
/S.uSco is the velocity-integrated CO line flux. This 
corresponds to a (frequency-integrated) luminosity of 
Lco(3-2) =3.2x10^ Lq, or a brightness-temperature lu- 
minosity of icO(3-2) ~ ^ ^^^^ ^ ^™ P*-^- 

The central frequency of the line is at -16 km s"^, cor- 
responding to a frequency shift of 4.9 MHz (so a central 
frequency of 91.801 GHz), so the redshift of the CO line 
is zco = 2.7668. At our signal to noise, the line shape 
is well described by a gaussian, there is no double pro- 
file visible, which could suggest a disk. There is also no 
hint of two separate lines which could suggest that the 
molecular gas is undergoing a major merger. 

In the case of AMS16, no line was detected (see Fig- 
ure 10). The full data cube was searched by looking 
for line emission at various spectral resolutions, using 
a routine provided by Roberto Neri, but no significant 
emission was detected at the position of the source. 

A signal at the 4.5 a level can be seen at 17 19 44 
+58 47 00. If real this would correspond to a line offset 
from the narrow-line redshift by ~600 km s"^ and with 
a full width at zero intensity of ~850 km . However, 
the position does not correspond to any source detectable 
in the R-band, 3.6 /im or 1.4 GHz image. This position 
corresponds to a 17 arcscc offset, or 116 kpc at the red- 
shift of the source. We do not consider this emission any 
further. 

The lack of detection of the CO (4-3) line can either 
be due to the line being shifted outside of the observed 
bandwidth (> ±1680 km s^^ shift), or the line being 
within the bandwidth but too faint. The redshift was 
determined from three narrow lines, Ly a (1216 A), N IV 
(1240 A) and C IV (1549 A) with individual redshifts of 
z = 4.171, 4.166 and 4.170 respectively (see Martmez- 
Sansigre et al., 2006a). Thus the maximum redshift offset 
of any individual line, with respect to the average {z) = 
4.169 is that of the N IV line with Az = 0.003, which 
only corresponds to a 175 km s"'^ offset. 

Given the bandwidth used, we consider the possibility 
that the CO (4-3) line falls outside of the observed range 
unlikely, so detection of CO (4-3) in AMS16 can be ruled 
out down to a 3a limit of Sco^'^ <0.23 Jy km s"^, 
assuming a box-car line shape with full width to zero 
intensity of 400 km s~^. From these limits, we now esti- 
mate the limits on the CO (4-3) luminosities and derived 
molecular gas properties. Figure 10 shows the spectrum 
extracted at the position of AMS16. 

At z = 4.169, ScoAi/ <230 mJy km s'^ cor- 
responds to a frequency-integrated line luminosity of 
Lco(4-3) <3xl0^ L0, or a brightness-temperature lu- 
minosity of ico(4-3) < 1 ^ 10^° K km s~i pc^. 

6.2. Inferred physical properties 

If we assume the molecular gas responsible for the 
CO emission is dense enough and warm enough for the 
(3-2) or (4-3) transitions to be thermalised, then the 

brightness-temperature luminosities of these high- J tran- 
sitions are the same as that of the (1-0) transition: 

■^CO(l-O) = -^CO(3-2) ^^'^ -^CO(l-O) = -^CO(4-3)- 

Following Downes & Solomon (1998), we can then esti- 
mate an upper limit on the mass of molecular hydrogen, 
by using 
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Ml 



H2 



Mr. 



L' 



= aco 



CO(l-O) 
Lq 



(10) 



where aco =0.8 (K km Mq)"-'^ is the conversion 
factor for local ULIRGs (Downes & Solomon 1998). 
Thus, under the assumption that ^00(1-0) ~ -^00(3-2)' 
wo infer a total gas mass of = 1-9 x 10"*^° Mq 

for AMS12. Assuming that i'co(i-o) ~ -^CO(4-3)' ^® 
can derive a limit for the total H2 mass in AMS16 of 

< 8 X 10" Mq. 

The far-infarcd hmiinosity of AMS12 of 2.5 x lO^^ 
(see Table 4), suggests a SFR of 4300 Mq yT~^. The ratio 
of far to CO luminosity is logio[LFiR/ico] — 3.0, the 
gas-to-dust ratio, Mg/Md, is 19, and the gas depletion 
timcscalc is estimated to be Tg ^4 Myr. 

For AMS16, we can only infer limits for those val- 
ues: if we assume the values from the T = 35 K, 
/? = 1.5 gray body, we infer logio[I'FiR/-f'co] ^2.4, 
Mg/Md<8 and Tg <16 Myr. Assuming instead the 
T = 47 K, /3 = 1.6 gray body yields the following values: 
logio[LFiR/ibo] >2.9, Mg/Md<20 and Tg <5 Myr. 

For both AMS12 and AMS16, the estimated values of 

logio[I'FiR/-f'co] "^-^^ consistent with 

the observed values for other higli-redshift sources with 
molecular gas measurements (e.g. Solomon & Vanden 
Bout 2005; Greve et al. 2005). The gas-to-dust ratios 
are much lower than the values ^100-150 found for local 
spiral galaxies. The values of Tg are ultimately derived 
from the quantity LprnJ L'^q and are thus also at the low 
end of the observed values. A possible explanation for 
the high values of Lfir/-^co that a fraction of the far- 
infrared luminosity is due to AGN heated dust, rather 
than by dust heated by young stars. 

We warn that if the assumption that the (4-3) tran- 
sition is thermalised is inappropriate, Li 



CO(l-O) 



> 



^co(3-2) or ^co(i-o) > -^co(4-3)' ^c wih have under- 
estimated Mg, and hence Tg and the gas-to-dust ratio. 
For a discussion about whether this assumption is ap- 
propriate, see Papadopoidos et al. (2001) and Ricchors 
et al. (2006), although for the high-redshift AGN stud- 
ied in detail it seems to be a reasonable assumption (see 
e.g. Beelen et al. 2004; Klamer et al. 2005; Aravena et al. 
2008). 

From the FWHM of the CO(3-2) hnc in AMS12, 
Al/ =275 km s~^, and assuming a characteristic radius, 
we can estimate the dynamical mass given (see Neri et al. 
2003): 



Mdyn sin^i , . , „4 
^— = 4 x 10 

Mq 



km s 



The CO is unresolved in the ~5 arcsecond resolution ob- 
servations with the PdBI observations, yielding a size 
limit of r <40 kpc. A better estimate for r is 2 kpc (Tac- 
coni et al. 2006). The inferred mass is then Mdynsin^i = 
6 X 10" Mq, and is obviously dependent on the unknown 
inclination angle i. The ratio of gas to dynamical mass 
is then Mn^/M^yn = 3 x sin^i. Thus the estimated gas 
mass is larger than the estimated dynamical mass, un- 
less the inclination angle i is < 35. It seems likely that 
the gas mass represents a high fraction of the total dy- 



namical mass, as found for other high-redshift galaxies 
observed in CO (Greve et al. 2005; Tacconi et al. 2006). 

7. COMPARISON TO OTHER SAMPLES OF Z ~ 2 

QUASARS 

Other samples of z ^ 2 quasars have been observed at 
millimetre or submillimetre wavelengths, allowing us to 
compare the far-infrared luminosities inferred for the dif- 
ferent samples. Far-infrared emission being isotropic, the 
properties of obscured and unobscured quasars should be 
similar if obscuration is only an orientation effect. To 
test this, we consider here three other studies at approx- 
imately the same redshift: 

• The sample of Omont et al. (2003, hereafter O03), 
consisting of optically-selected unobscured quasars 
observed at 1.2 mm. These sources arc intrinsically 
very bright with Mb ^ —27, corresponding to Lboi 
in the range 5 x 10^^ - 5 x 10^^ Lq, and at 1.94 < 
z < 2.85. They have a detection rate of 9 out of 26 
sources (35%). 

• The sample presented by Page et al. (2004, P04), 
which consists of optically unobscured quasars 
known also to be X-ray unabsorbed. These sources 
have values of Ly^oi^ 10^'^ Lq, similar to our sam- 
ple, with 1.20 < z < 2.35. They were observed at 
850 /um, with a detection rate of 1 out of 20 (5%). 

• The optically unobscured but X-ray absorbed 
quasars from the sample of Stevens et al. (2005, 

S05), with similar values of Lboi to our sample, 
1.01 < z < 2.80. They were observed at 850 fim 
with 6 out of 19 sources detected at the > 3a level 
(32%). 

In all samples, the bolometric Imninosities have been 
estimated consistently using the Elvis et al. (1994) SED, 
and correcting for absorption at X-ray energies (S05) 
or dust-extinction at optical wavelengths (this sample). 
The far-infrared luminosities have all been computed as- 
suming a gray body with T = 47 K, /? = 1.6. Changing 
these parameters will affect very similarly the far-infrared 
luminosities of all samples. 

Figure 11 shows ^fir vs Lboi for the four samples of 
2; ~ 2 quasars. The coloured symbols represent the indi- 
vidual sources, the black filled symbols the mean values 
for each sample, and the black empty symbols the mean 
of the non-detections. The dashed line shows a constant 
far-infrared fraction of ^^^^ = 0.05, characteristic of the 

low-redshift quasars studied by Elvis et al. (1994). 

From the available data, the O03 and P04 samples as 
a whole have a ratio of similar to the low-redshift 

-t>bol 

quasars. However, our sample and the S05 samples both 
seem to have larger far-infrared fractions than the O03 
and P04 samples^^. 

It is difRcult to quantify accurately the statistical significance 
of the differences between samples. As quoted in the ASURV man- 
ual: "It is possible that all existing survival methods will be inac- 
curate for astronomical data sets containing many points very close 
to the detection limit." (Isobe & Feigelson 1990), which is clearly 
the case in this comparison. To make things worse, the individual 
samples considered here are small in a statistical sense: they all 
consist of less than 30 objects. 
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If all obscuration were orientation dependent then the 
host galaxies of obscured and unobscured quasars would 
be expected to have, on average, the same properties, 
such as star-formation rates (e.g. Clavel et al. 2000) and 
cool dust content (e.g. Haas et al. 2004). If, however, 
some obscured quasars are at a different evolutionary 
stage, then their host galaxy can be reasonably expected 
to have a larger gas and dust content, which might be 
causing the obscuration, or at least contributing to it. 
Obscured quasars might therefore be in a different evo- 
lutionary phase to unobscured quasars (e.g. Sanders et al. 
1988; Fabian 1999). 

The extinction due to dust is negligible in the far- 
infrared (FIR), so that at wavelengths ^ 50 ^m the emis- 
sion becomes essentially isotropic. Hence, if orientation 
is the only difference between obscured and unobscured 
quasars, they should have virtually identical cool dust 
properties, irrespective of whether this emission is due 
to dust heated by the quasar or by young stars in the 
host galaxy. 

Although it is difficult to quantify the significance of 
the difference between samples, it seems that optically 
obscured quasars (and optically-unobscured but X-ray 
absorbed quasars) have, on average, higher far-infrared 
luminosities than unobscured quasars. In turn, this im- 
plies that the cool-dust content of the host galaxies are 
different so that, on average, obscured quasars are hosted 
by dustier galaxies. 

This is, of course, partly a result of selection. The sam- 
ples selected to contain obscured sources will of course be 
biased towards dustier objects, but such a bias will only 
show in the far-infrared observations if obscuring dust 
on kpc scales is important (as opposed to dust in the 
torus only, which will emit mostly at mid-infrared wave- 
lengths). The S05 sample is particularly biased towards 
such host obscuration, since it contains optically-bright 
quasars that are not obscured by the torus, but which 
show absorption at X-ray energies. 

Hence, the difference in far-infrared luminosities be- 
tween obscured and unobscured quasars lends strong sup- 
port to the scenario that some of the obscuring dust is 
in the host galaxy, not only the torus, and there is a hint 
that many obscured quasars might represent a different 
evolutionary phase to unobscured quasars. 

8. DISCUSSION AND CONCLUSIONS 

We have observed a sample of z 2 radio-intermediate 
obscured quasars at 1.2 mm (250 GHz) using MAMBO. 
The detection rate is 5 out of 21 (24%). Sources with 
narrow lines seem to have a higher detection rate than 
sources with no narrow lines, but with such small num- 
bers the difference is not significant. Larger samples will 
be required to study any differences. 

Stacking leads to a statistical detection of {Si,2 mm) = 
0.96 ± 0.11. Even if only the non-detections are stacked, 
they still yield a statistical detection, with (5i.2 mm) = 
0.51 ± 0.13. Thus, the typical flux density of this sample 
is ^0.5-1.0 mjy, and this corresponds to a far-infrared 
luminosity ^ 4 x 10^^ L©. 

If the far-infrared luminosity is powered entirely by 
star-formation, and not by AGN-heated dust, then the 
typical inferred star-formation rate is ~700 M© yr"-'^. 
This large star-formation rate is comparable to those in- 
ferred in ULIRGs and SMGs, the most powerful star burst 



galaxies known. 

The observations at 1.2 mm also allow us to estimate 
the mass of the cool dust, and we find a typical mass 
of ~ 3 X 10® Mq. If heated by young stars, this dust 
is expected to be distributed on ~2 kpc scales (the scale 
found for SMGs). If it is heated by the central AGN, with 
Lhoi^ 10^^ Lq, then it is expected to be distributed on 
scales ^10 kpc. 

Indeed, we estimate that such a large mass of cool dust 
is capable of causing alone the large values of Ay in this 
sample, without help of an obscuring torus. 

Combining our observations at mid-infrared and mil- 
limetre wavelengths, we present dust SEDs for our sam- 
ple, and derive a typical SED for our sample of high- 
redshift obscured quasars. The clumpy tori models can- 
not reproduce the SED. However, clumpy tori models 
with an additional screen of cold dust (and the corre- 
sponding far-infrared emission) do reproduce this mean 
SED. The amount of extinction from this additional 
screen is derived from the cool dust mass, from the ob- 
servations at 1.2 mm. The depth of the silicate feature 
can be consistently achieved by the inferred cool dust 
mass provided that the dust is on scales ^2 kpc, which 
is in excellent agreement with the typical radius of the 
far-infrared emission in SMGs. Again, this lends sup- 
port to kpc-scale dust along the host galaxy playing an 
important role in the obscuration of these sources. 

Obscuration by dust on kpc scales would also explain 
why in about half of the sample the narrow emission lines 
from the central AGN are not detectable, as well as why 
in some cases the jets seem to be pointing towards us, 
something unexpected if they are only obscured by the 
torus of the unified schemes. 

However, we remind the reader that unobscured 
quasars have also been found to have similar cool-dust 
masses, and that in our sample there is not a one-to-one 
correspondence between the presence (or lack of) narrow 
lines and detection at 1.2 mm. The presence of a large 
mass of dust does not guarantee obsciiration. A clumpy 
dusty interstellar medium where different lines of sight 
have vastly different optical depths is a likely explana- 
tion for why some quasars with large dust masses are 
obscured while others are not. 

When comparing to other samples of high-redshift 
quasars, we find that the obscured quasars probably have 
a higher fraction of their luminosity emerging at far- 
infrared wavelengths, compared to unobscured quasars. 
The far-infrared emission is isotropic, so that this dif- 
ference cannot be ascribed to orientation-dependent ob- 
scuration: obscured quasars are likely to have higher 
far-infrared luminosities and cool-dust masses than un- 
obscured quasars, suggesting the host galaxies are in a 
dustier, presumably earlier, phase. 

Additionally, we have searched for molecular gas in 
two sources a.t z = 2.767 and 4.169, using the CO (3-2) 
and (4-3) transitions. In the z = 2.767 source we de- 
tect a line with Lco(3-2) =3.2x10^ Lq (equivalent to a 
brightness-temperature luminosity of -^co(3-2) ~ '^■^ ^ 
10^° K km s~^ pc^). In the other source, the lack of de- 
tection suggests a line luminosity Lco(4-3) <3xl0^ Lq 
(-^^00(4-3) < 1 X 10^° K km s~^ pc^). Under the as- 
sumption that in these objects the (3-2) and (4-3) tran- 
sitions are thermalised, we can estimate the molecular 
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gas contents to be Mr^ = 1.9 x 10^" and < 8 x 10^ Mq, 
respectively. The estimated gas depletion timescales are 
Tg = 4 and <16 Myr, and low gas-to-dust mass ratios of 
Mg/Md= 19 and < 20 are inferred. A dynamical mass of 
Afdynsin^i = 6 X 10^ Mq is estimated from the CO(3-2) 
detection. 
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Table 1 

Observed properties of the sample of obscured quasars. 



Name 


RA 
(J2000) 


Dec 




''■spec 


Optical'' 


Radio'^ 


■/core 


Mid-infrared*^ 


Si. 2 mm'' 

/mJy 


AMSOl 


17 


13 


11.17 


+59 


55 


51, 


,5 


_ 


B 


SS 


<0.3 


C 


0.63+0.53 


AMS02 


17 


13 


15.88 


+60 


02 


34, 


,2 


1.8 


B 


ss 


_ 


S 


0.87+0.54 


AMS03 


17 


13 


40.19 


+59 


27 


45, 


,8 


2.698 


NL 


SS 


0.8 




0.76+0.51 


AMS04 


17 


13 


40.62 


+59 


49 


17, 


,1 


1.782 


NL 


ss 


- 


• 


0.55+0.58 


AMS05 


17 


13 


42.77 


+59 


39 


20, 


,2 


2.850 


NL 


GP 


1.1 


S 


-0.26+0.40 


AMS06 


17 


13 


43.91 


+59 


57 


14, 


,6 


1.8 


B 


ss 


<0.3 


s 


0.86+0.41 


AMS07 


17 


14 


02.25 


+59 


48 


28. 


,8 


_ 


B 


FS 


_ 


c 


1.05+0.44 


AMS08 


17 


14 


29.67 


+59 


32 


33. 


,5 


1.979 


NL 


- 


- 


s 


1.99+0.34 


AMS09 


17 


14 


34.87 


+58 


56 


46. 


,4 


2.1 


B 


ss 


0.9 


s 


0.51+0.56 


AMSIO 


17 


16 


20.08 


+59 


40 


26. 


,5 


- 


B 


ss 


- 


c 


0.62+0.79 


AMSll 


17 


18 


21.33 


+59 


40 


27, 


,1 


1.6 


B 


ss 




P,S 


0.36+0.57 


AMS12 


17 


18 


22.65 


+59 


01 


54, 


,3 


2.767 


NL 


ss 


0.3 




3.74+0.60 


AMS13 


17 


18 


44.40 


+59 


20 


00, 


,8 


1.974 


NL 


ss 




s 


0.27+0.519 


AMS14 


17 


18 


45.47 


+58 


51 


22, 


,5 


1.794 


NL 


ss 




s 


0.27+0.44 


AMS15 


17 


18 


56.93 


+59 


03 


25, 


,0 


2.1 


B 


FS 


0.5 


s 


0.12+0.70 


AMS16 


17 


19 


42.07 


+58 


47 


08, 


,9 


4.169 


NL 


GP 


<0.6 




2.12+0.54 


AMS17 


17 


20 


45.17 


+58 


52 


21, 


,3 


3.137 


NL 


SS 


<0.3 


P,S 


1.73+0.48'' 


AMS18 


17 


20 


46.32 


+60 


02 


29, 


,6 


1.6 


B 


FS 




s 


1.11+0.53 


AMS19 


17 


20 


48.00 


+59 


43 


20, 


,7 


2.3 


B 


FS 


1.0 


s 


2.76+0.66 


AMS20 


17 


20 


59.10 


+59 


17 


50. 


,5 




B 


GP 




c 


-0.78+0.93 


AMS21 


17 


21 


20.09 


+59 


03 


48. 


,6 


1.8 


B 


SS 


0.7 


P,S 


0.52+0.65 



''Rcdshifts with 3 decimal places are from optical spectroscopy (Martfnez-Sansigre et al. 2006b), while those with only 1 decimal place 
are from mid-infrared spectroscopy (Marti'ncz-Sansigre et al. 2008). The redsliift for AMS05 is from the optical spectrum taken by Smith 
et al. (2009). 

''Summary of optical spectroscopy properties. B stands for blank spectrum, NL for narrow lines. 

"Summary of radio spectral properties from Martmez-Sansigre et al. (2006a). SS stands for steep spectrum, FS for flat spectrum and GP 

for gigahertz-peaked source. 

''Approximate fraction of radio flux at 1.4 GHz recovered at VLBI resolution at 1.6 GHz (Klockner et al. 2009). 

"Spectral properties in the mid-infrared (see Martmez-Sansigre et al. 2008). C stands for continuum only, S for silicate absorption, P for 
polycyclic aromatic hydrocarbons (PAHs). 
'Flux densities at 1.2 mm. Significant detections are marked in bold, 
spiux density for AMS13 taken from Lutz et al. (2005), their source 22204. 
•"Flux density for AMS17 taken from Sajina et al. (2008), their source 22558. 



Table 2 

Comparison of resulting fluxes in Win07 and Sum08/Win08. Bold shows the values used in Table 1 



Name 


Win07 


Sum08/Win08 


Combined 




/mJy 


/mJy 


/mJy 


AMS07 


2.43+0.83 


1.05+0.44 


1.27+0.41 


AMS08 


2.04+0.59 


1.99+0.41 


1.99+0.34 


AMSIO 


0.81+0.93 


0.45+1.1 


0.62+0.79 


AMSll 


0.83+0.64 


-0.59+1.00 


0.36+0.57 


AMS12 


4.03+1.03 


3.74+0.60 


3.80+0.55 


AMSl?" 


2.39+0.55 


1.45+0.62 


1.94+0.44 


AMS20 


-0.16+1.08 


-1.24+1.56 


-0.78+0.93 



^The flux density for AMS17 used in Table 1, 1.73+0.48, is taken from Sajina et al. (2008), their source 22558. 
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Table 3 

Inferred properties of the sample of obscured quasars. 



Name logiofLbol/L©]" ^v" T = 47 K jS = l.e*" T = 35 K /3 = 1.5" 
logiofLFiR/Lp] logio[SFR/M0 yr-^] logio[Md/MQ] logiolLpiR/LQ] logio[SFR/MQ yr"!] logio[Md/MQ] 



AMSOl 


>13.2 


> 40.5 


<13.0 


< 3.2 


< 8.7 


<12.4 


< 2.6 


< 8.9 


AMS02 


12.4 


0.5 


<13.0 


< 3.2 


< 8.7 


<12.4 


< 2.6 


< 8.9 


AMS03 


>13.2 


> 21.5 


<12.9 


< 3.2 


< 8.6 


<12.3 


< 2.6 


< 8.9 


AMS04 


13.7 


68.5 


<13.0 


< 3.3 


< 8.7 


<12.4 


< 2.7 


< 9.0 


AMS05 


13.6 


5.5<* 


<12.8 


< 3.0 


< 8.5 


<12.2 


< 2.5 


< 8.8 


AMS06 


13.1 


12.0 


<12.9 


< 3.1 


< 8.6 


<12.3 


< 2.5 


< 8.8 


AMS07 


>12.9 


> 17.0 


<12.9 


< 3.1 


< 8.6 


<12.3 


< 2.5 


< 8.8 


AMS08 


>13.6 


> 65.0 


13.1 


3.3 


8.8 


12.5 


2.7 


9.0 


AMS09 


13.1 


16.0 


<13.0 


< 3.2 


< 8.7 


<12.4 


< 2.6 


< 8.9 


AMSIO 


>13.1 


> 55.0 


<13.2 


< 3.4 


< 8.9 


<12.6 


< 2.8 


< 9.1 


AMSll 


12.7 


13.0 


<13.0 


< 3.3 


< 8.7 


<12.4 


< 2.7 


< 9.0 


AMS12 


13.3 


21.0 


13.3 


3.5 


9.0 


12.7 


3.0 


9.3 


AMS13 


13.9 


22.5 


<13.0 


< 3.2 


< 8.7 


<12.4 


< 2.6 


< 8.9 


AMS14 


>13.8 


> 73.5 


<12.9 


< 3.2 


< 8.6 


<12.3 


< 2.6 


< 8.9 


AMS15 


>13.0 


> 40.0 


<13.1 


< 3.3 


< 8.8 


<12.5 


< 2.8 


< 9.1 


AMS16 


13.9 


21.0 


12.9 


3.2 


8.6 


12.4 


2.7 


9.0 


AMS17 


13.8 


33.5 


12.9 


3.2 


8.6 


12.4 


2.6 


8.9 


AMS18 


>14.2 


> 100.0 


<13.0 


< 3.2 


< 8.7 


<12.4 


< 2.6 


< 8.9 


AMS19 


13.4 


7.5 


13.2 


3.4 


8.9 


12.6 


2.9 


9.2 


AMS20 


12.7 


0.2 


<13.2 


< 3.5 


< 8.9 


<12.6 


< 2.9 


< 9.2 


AMS21 


13.1 


24.0 


<13.1 


< 3.3 


< 8.8 


<12.5 


< 2.7 


< 9.0 



''Bolometric luminosities and extinctions derived from the broad-band data between 3.6 and 24 /im, as described in Section 4. 

'^Far-infrared luminosities, star-formation rates and cool-dust masses inferred from the 1.2 mm flux densities, assuming a gray body with 
T = 47 K and /3 = 1.6 (sec Section 3). 

■^Far-infrared luminosities, star-formation rates and cool-dust masses inferred from the 1.2 mm flux densities, assuming a gray body with 
T = 35 K and /3 = 1.5 (see Section 3). 

■^The value of Av =5.5 was inferred by Smith et al. (2009) using additional data at 1.2, 1.7 and 2.2 //m. 



Table 4 
Summary of gray body fits 



Source 


13 T 


loglO [ipiR 


logio[Md 




/K /Lq] 


/M0] 




AMS08'' 


1.1 57 


13.0 


8.9 


AMS12'' 


1.5 52 


13.4 


9.0 


AMS19'' 


1.5 42 


12.9 


9.0 



"In the case of AMS08, three data points are used: 70 and 160 nm as well as 1.2 mm, so that three parameters can be fitted: /?, T and 

^FoT these two sources, the 70 fj,ui data do not provide useful constraints, and only two parameters are fitted: T and ipiRi while /? is 
kept fixed at a value of 1.5. 
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Fig. 1. - Illustration of the physical components that dominate at each observed band. The figure shows the SED of a model torus 
(solid line, from Ncnkova ct al. 2008a, b) and a cool dust component (dashed line) modelled as a gray body with temperature T =40 K 
and cmissivity index /3 = 1.5 (sec Equation 2). The observed bands at 24, 70 and 160 ^m, as well as 1.2 mm, are available for all 21 
sources and have been marked in gray. This model source is assumed to be at z = 2 so the bands have been accordingly shifted to the 
rest-frame frequencies. The relative luminosities of the two components arc arbitrary but representative of the sources in this sample. For 
completeness, the exact torus parameters used arc; single-cloud optical depths tv = 60, gaussian angular distribution of clumps around 
the equator, a = 60, the number of clouds along a ray in the equatorial plane N = 10, the ratio of inner to outer torus radius Y = 30, and 
the power of the radial density distribution, g = 0. 
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Fig. 2. — Flux density, Si. 2 mm, versus redshift 2 for the sample of obscured quasars. Sources marked in blue have redshifts from optical 
spectroscopy, those marked in red have redshifts only from their mid-infrared spectra. All limits are 3(t. The four green sources have no 
spectroscopic redshifts and have been placed at 2 = 1.98, the median value for sources with a spectroscopic redshift. At 2 = 2 there are 4 
green sources (one detection and three upper limits) and 2 blue sources (one detection and one upper limit). The dotted line represents 
the approximate flux density limit of the MAMBO observations: ~1.65 mjy. The dashed line shows the value of (Si. 2 mm), the mean flux 
density for all 21 sources. 
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Fig. 3. — Distribution of observed flux densities. The empty histogram represents the measured flux densities for all sources, the shaded 
histogram for detections. The distribution is clearly centred on a positive value, and skewed towards the high-flux density end. The solid 
vertical line shows the stacked mean of all sources, the dashed line the stacked mean for detections, and the dotted line the stacked mean 
for non-detections. 
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Fig. 4. — Far-infrared luminosity, Lpia, versus redshift z, for the two fiducial gray bodies. Symbols are the same as in Figure 2. The 
dotted line represents the value of LpiR corresponding to the approximate flux limit, while the horizontal dashed line represents the 
far-infrared luminosity corresponding to (5i.2 mm) at each redshift. 
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Fig. 5. — Spectral energy distributions of the sample of high-redshift obscured quasars. Data points are marked as black points for 
detections or arrows for limits (left to right): 1.2 mm, 160, 70, 24, 8, 5.8, 4.5, 3.6 fim (data from Lacy et al. 2005; Fadda et al. 2006; 
Frayer et al. 2006, and deeper unpublished IRAC data). All limits are 3cr. Overlayed in black are also the IRS spectra, when available. 
The far-infrared BED is modelled by a gray body, as described in Section 3.2. The mid-infrared component is modelled by an unobscured 
quasar SED with foreground extinction (see Section 4). The dashed lines are the Elvis et al. (1994) quasar SED with a screen of dust (with 
the value of Ay as quoted in Table 3 and using the MW dust model of Pei 1992) and a, z = elliptical galaxy from Coleman et al. (1980). 
The quasar SED has been truncated at wavelengths longer than rest-frame 80 fim. The blue solid line is the sum of both obscured quasar 
and host galaxy SEDs at wavelengths shorter than 20 fim, fitted using the data points between 3.6 and 24 fim (see Section 4). Source 
AMS05 has additional data points at 1.2, 1.6 and 2.2 fim from Smith et al. (2009). 
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Fig. 5. — Continued. Source AMSll has an extra point at 850 fim, from the observations of Frayer et al. (2004), while source AMS12 
has a hmit on the 3 mm flux density from the PdBI observations. Sources AMS08 and AMS12 have an extra gray body overplotted (black 
solid lines), described by the values in Table 4. 




Fig. 5. — Continued. Sources AMS13 and AMS16 have additional data at 350 fim from the SHARC observations, 
has a limit on the 3 mm flux density from the PdBI observations. 



In addition, AMS16 




Fig. 5. — Continued. Source AMS19 has an extra gray body overplotted (black solid line), described by the values in Table 4. 
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1012 IQ13 1014 

Fig. 6. — Mean dust SED for the sources in our sample with spectroscopic redshifts (solid red lines, see Section 5). The near-infrared 
component is the mean of the observed IRAC photometry (u = 1 — 4 X 10^'' Hz, observed at 3.6-8.0 /^m), the mid-infrared component 
(y = 2 — 6 X 10^'' Hz) is constructed using the IRS spectra, while the far-infrared component {u < 1 X 10^^ Hz) described by either one of 
the two fiducial gray bodies (the top solid red line has T =47 K, (3 =1.6, the bottom one has T =35 K, /3 =1.5). Overplotted for reference 
are an unobscured quasar SED (dotted line Elvis et al. 1994) as well as the same quasar SED with a screen of MW dust (from Pei 1992) 
with the mean Ay= 32 (plotted as a dashed line), and a 2 = elliptical galaxy Coleman et al. (dash-dotted line, from 1980). This last 
SED has been normalised to be the progenitor of a present day ^ 2L* galaxy. 
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Fig. 7. — Comparison of the mean SED (solid black lines) with different clumpy models of AGN tori. The progenitor of a 2L* galaxy has 
been overplotted again (dot-dashed line). The clumpy models, from Nenkova et al. (2008a, b), all have single-cloud optical depths Ty = 60. 
Left panel: Different lines show models with varying parameters for the width of the gaussian angular distribution of clumps around the 
equator, a, the number of clouds along a ray in the equatorial plane N, the ratio of inner to outer torus radius Y, and the power of the 
radial density distribution, q. The models cannot simultaneously reproduce the relatively flat continuum between 1 X lO^^'' and 7 X lO^^'' Hz, 
the sudden drop at higher frequencies, the depth of the silicate feature and the millimetre detection. Right panel; Additional extinction 
from dust along the line of sight has been added to the two flatter clumpy models from the left panel. The emission from this cool dust is 
represented by the two fiducial gray bodies used throughout this work. The typical quasar SED of our sample can be reasonably reproduced 
using the clumpy torus models, with foreground extinction Ay= 80 and re-emission from a mass of 3 X 10 Mq of cool dust. The emission 
from this cool-dust mass and the foreground extinction are consistent provided the dust is within a 2 kpc radius (using Equation 8 from 
Section 3.6). 
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Fig. 8. — Comparison of the mid-infrared component of the mean SED (black sohd hne) with four local ULIRGs: IRAS 12514+1027 
(hereafter IRAS 12514), Mrk 231, Mrk 273 and Arp 220. 
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Fig. 10. — Spectrum of AMS16 obtained with the Plateau de Bure Interferometer. The frequency was centred on 89.2002 GHz, the 
expected location of the CO (4-3) rotational transition. The rms noise per 10 MHz bin (33 km s"!) is 0.65 mjy beam"!, assuming a 
box-car CO line shape with full width to zero intensity of 400 km s~^ and the line to be within the frequency covered, the 3ct limit is 
Ai^Sco <230 mJy beam"! km s'^, or ^^0(4-3) < 1 ^ 10^° ^ km s'^ pc^ 
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Fig. 11. — ipiR vs Lbol for the four samples of z ^ 2 quasars considered. O03 is the sample presented by Omont et al. (2003), P04 
by Page et al. (2004) and 305 by Stevens et al. (2005). Black filled symbols represent the weighted mean of all sources in each sample, 
black empty symbols represent the weighted mean of non-detections in each sample. The dashed line represents Lfir= 0.05Lboli the value 
found in low-redshift quasars considered by Elvis et al. (1994). For reference, the line LpiR=Lbol has been plotted as a dotted line. For 
all samples, a gray body with T = 47 K and /3 = 1.6 has been assumed to calculate LpiR. 



